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Abstract 


The absolute Intensities, positions, rotational constants and half-widths 

-1 

of the lines of three overlapping bands of CO near 4640 cm have been 

2 

determined by nonlinear, least-squares, fitting procedures. The intensities 

0 0 1 1 16 12 18 
of (00 2) «- 00 0 and (01 2) «- (01 0) of 0 C 0 are In good agreement 

I II 

with the estimates in the AFGL listing. An improved value has been obtained 

1 0 16 12 16 

for the intensity of (31 0) «- 00 0 of 0 C 0. 

Ill 
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1. INTRODUCTION 

Carbon dioxide is one of the more important, atmospheric, infrared- 

absorbing gases due to its relatively high and increasing concentration. The 

spectral parameters of its bands are required for understanding radiative 

transfer in the atmosphere and for many other applications. In this paper, 

the line intensities and self-broadened widths as well as the band centers of 

three overlapping bands of CO near 2.2 pm are reported. These were obtained 

2 

by nonlinear, least-squares (NLLS) regressions based on the method of Lin et 
al . ( 1 ) as modified by Hoke and Shaw (2,3). 

2. DATA 

The experimental conditions of the CO samples are summarized in Table 

2 

-1 

I. The bands occur between about 4545 and 4662 cm , Analyses of similar 

-1 

bands near 2600 cm in these sajne spectra were reported recently (2,3) and 

these references give descriptions of the spectra and the calibration 

procedures. Because of the higher detectivity of the InSb detector at longer 

wavelengths the SNR of the spectra previously analyzed was in excess of 500. 

Fig. 1 shows a portion of Spectrum I. The main band (centered around 
-1 

4640 cm ) has no Q-branch. Because there is a weak unresolved Q-branch near 

-1 -1 
4591 cm , the region between 4590 and 4592 cm was assigned a weight of zero 

in all of the analyses. The data consisted of signal values obtained from 


unapodized interf erograms . 
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3. THEORY 

A spectrum similar to the observed spectrum is calculated from model 
functions with spectral quantities as adjustable parameters. In this analysis 
the models consisted of the upper and lower state energy levels, the line 
strengths and the dependence of the collision-broadened half-widths on 
rotational quantum number. The influences of the source and the observing 
instrument on the spectrum were also taken into consideration by modeling the 
effects of the retardation A of the interferometer mirror, the uncorrected 
phase error $ and the angular spread of the rays passing through the optical 
system. The frequency dependence of the baseline was described by a single 
parameter but, for the background, a five-degree polynomial was needed. 

The calculated signal values S*(v) are convolved with an instrument 
spectral response function ISRF(v) as expressed by the relation: 

’(v) ISRF(v^,v) dv. [1] 


S( V - I s 


The spectral response function used in this analysis is of the form (4), 


ISRF ( A , 4> , v , d ) = ¥<d) 


Acos«(> sinc(2uA) + sin4> sin (uA)/uA 


[ 2 ] 


where u = v(v - v ), and the sine and cosine functions model the uncorrected 
o 

phase error <|>. ’i'(d) is a rectangular function whose half-width d depends on 

the angular spread of the rays passing through the interferometer. This has 
been described as the ’extended source effect’ (4). Unlike the retardation, 
the phase error and extended source parameters are frequency dependent. 

The line positions are given by the difference of term values 


v 

o 


(cm 1 


F’(v’,J’,l’) - F"(v",J",l"), 


C3] 
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where the primes and double primes Indicate the upper and lower levels 
respectively. The term values are 


F(v,J,l) = Q + B yJl [J(J+l) - l 2 ] - D vJl CJ(J+l) - l 2 ] 2 + , [4] 

where the G's are the vibrational terms and the B*s and D’s are the rotational 
parameters. 

Individual line intensities are expressed (5, 6), by; 

S(m) = [8 v 3 C |R(0) | 2 /3hcQ_] v (m)F(m)C (m) exp[-E"(m)/kT] , [5] 

N It O O 

where m = -J for the P-branch and m = J + 1 for the R=branch, and 


2 

F(m) = (1 + f^m + B ? m + ... ) 
C Q (m) = (m 2 - l 2 )/ |m| 

F(m) = (1 + pa) 2 
C Q (m) = |m + 1 |/2 


for the * - « and 
£ - X bands. 

for the w - £ band. 


C includes the vibrational partition sum and the natural abundance of the 
N 

isotope. The other symbols have their usual meaning. This relationship 

yields intensities which can be compared directly with those in the AFGL 

Listing (7). A natural CO sample was used and the isotopic abundances are 

2 

assumed known. 

As expected, the Lorentz and Vbigt profiles yielded essentially identical 
results for the line widths and intensities. The |m| dependence of the widths 
was modeled by using the empirical relationship of Winters et al.(8). 


a(m) = <* 0 <P.T) 


V Q + W x exp(V 2 |m|) + W 3 |m| exp[-B"|m| ( |m| - l)/kT]J , 


[ 6 ] 


where o q (P,T) = P/P o (T q /T) , P Q = 760 torr, T q = 296 K, and the W^ are 
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regression parameters. This expression was found to be appropriate In the 
previous analyses (2). 

The regression program compares the calculated signal values with the 
experimental values. The set of parameters which gives the minimum value for 
the sum of squares of the differences between these signal values Is 
considered to give the best estimates for the spectral quantities. The total 
number of data points included in the analysis is 3883 and the number of 
parameters usually retrieved is between 30 and 40. 


4. RESULTS 

It is convenient to define a "slgnal-to-difference" ratio (SDR) by 
analogy with the slgnal-to-nolse ratio (SNR), where the "difference" Is 
defined as the standard deviation of the differences between the calculated 
and the experimental signal values. The best values of the SDR were 110 for 
Spectrum I and 220 for Spectrum II. The corresponding values of the SNR are 
130 and 260. The differences between the SDR and SNR are primarily due to 
very weak lines of unmodeled hot bands. 

The experimental and the calculated signal values of a portion of band I 
from Spectrum II are compared in Fig. 2. A plot of the residual is also shown 
below on an expanded scale. A ’noise spectrum’ obtained from a region with no 
significant absorption is also given together with its corresponding residual 
spectrum which, in this case, is the difference from a straight line through 
the noise. The appearances of these two residual spectra are very similar. 
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a. Vibration-Rotation Constants 

The vibration-rotation parameters in Table II were only determined from 
Spectrum I since the broader lines of Spectrum II are unsuitable for precise 
line-position-parameter determination. The numbers of lines modeled in each 
band are given in Table III. The quoted uncertainty in each parameter is 3 a, 
where a is the standard deviation estimated by the regression program. Two 
sets of rotational constants were obtained for band I and band II; one from a 
retrieval with no constraints and the other with the lower state rotational 
constants fixed to some recently reported values (7, 11, 12). In neither case 
was it possible to obtain significant estimates of the H's and these 
parameters have ‘been assigned the value zero. 

The band centers in Table II include the contribution due to the 
vibrational angular momentum quantum number, so that 

v „ = AG , „ + <B"!" 2 + D"4" 4 - B*4' 2 - D'i' 4 ) , [7] 

center v* ,v” 

where AG is the difference between the upper and lower vibrational 

v' ,v" 

levels. Other estimates of band parameters are also included in Table II. 

Different estimates for the relative line positions for band I are 

compared in Fig. 3. The differences are given by v - v where the v are the 

a r r 

line positions calculated with the constants obtained from the best 

unconstrained retrieval and v are the positions obtained by assuming a common 

a 

band center and either the parameters reported by other investigators or 
obtained from other retrievals assuming a common band center. The three 
unlabeled curves near the zero line represent parameters from three different 
unconstrained retrievals with essentially the same SDR values. Their 
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agreement is a measure of the internal consistency of the results. Line 
position differences for the other bands are less consistent since the bands 
are weaker and hence their parameter estimates are less precise than those of 
band I. 

The rotational parameters from Spectrum I were used to calculate the line 

positions of Spectrum II but new values for the band centers of Spectrum II 

-1 

were retrieved. These showed a shift of about 0.01 cm from the low pressure 
positions. Shifts of this magnitude are within the range routinely observed 
on spectra measured on different days with our equipment. They are attributed 
to small changes in the optical alignment of the spectrometer (2) and they 
make it necessary to calibrate each spectrum individually (13). 


b. 


Intensities 

The band intensity is the sum of the line intensities S^(v): 


S ° = l S <v) , 
i=l 


The dipole moment matrix elements and band intensities obtained from both 

spectra are given in Table IV. A strong Coriolis interaction in band II 

causes the Intensity asymmetry of the P- and R-branches shown in Fig. 4. Line 

Intensities for this band vanish near m = -6. The existence of similar 

asymmetries in other v - £ bands of CO has also been reported (14, 15, 16, 

2 

17). The subbands of band III were treated as separate bands with a common 
band center and dipole moment. The reported band intensity is the sum of the 
intensities of both sets of lines. It was found that vibration-rotation 
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interaction parameters could only be retrieved for band II. 

c. Half -widths 

The Lorentz widths as functions of |m| shown in Fig. 5 were obtained from 
the parameter values of Table V. There is good agreement for |m| > 10. The 
same width parameters were used for all bands in each spectrum. In Fig. 6 the 

widths from the high pressure spectrum are compared with values obtained from 

-1 -1 
the 2600 cm region of the same spectrum and from near 2760 cm in another 

spectrum recorded at 741 torr and 292 K (2, 18). There is some evidence of a 

weak dependence of the widths on the vibrational transition although the 

differences between the previously estimated widths and those reported in this 

work are within 3 a, where a is the standard deviation in a(m) estimated from 

the standard deviations of the parameters of Bq. (6). 

d. On background, baseline and Instrument effects 

A five-degree polynomial was required to model the background. The 
baseline was more difficult to determine because of the small absorptances of 
the lines (average peak absorptance value ~0.25). For Spectrum II a value of 
-0.100(2) was found. This is about 10% of the observed signal for zero 
absorptance. In Spectrum I the retrieved value of the baseline shift was 
smaller than its uncertainty and hence it was set to zero. The retrieved 
values for the retardation shown in Table VI agree with those estimated 


previously (4). 
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5. DISCUSSION 

Despite the decreased SNR the overall precision obtained in the 

determination of Intensities and rotational constants is comparable with that 

-1 

reported for the 2600 - 2800 cm region in the previous analyses. As 

f 

expected, the most precise constants and intensities were obtained from the 
strongest band (band I). The intensity of band II appears to be significantly 
greater than the estimate given in the AFGL listing (7). We are not aware of 
any other reported experimental values for these band intensities. 

The dipole moment matrix element values in Table IV agree within 1% for 
band I and 2 % for band II and band III . This is within the estimated 
experimental errors. 

It appears, from Table V, that, although there is some uncertainty in the 

width parameters, the CO widths from different spectral regions and different 

2 

spectra are essentially consistent within the experimental uncertainty. 
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Table I 


Experimental conditions of the CO samples analyzed. 


Condition 

Spectrum I 

Spectrum II 

Pressure (torr) 

177 

658.6 

Temperature (X) 

294.2 

294.2 

Pathlength (m) 

476 

131 

SNR 

-130 

-260 
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Table II 


Retrieved rotational constants compared with other reported values. The 
letters at the end of rows Identify the constants according to the following 
assignments; M: Mandin (Ref. 9), C: Chedln (Ref. 10), R: AFGL listing 

(Ref. 7), F: Lower state rotational constants fixed to the values reported by 

Guelachvlll (Ref. 11) for band I, by Sakai (Ref. 12) for band II and by AFGL 
(Ref. 7) for band III. The upper state constants are retrieved, U: Constants 
retrieved with no constraints. 

The figures In parentheses are the uncertainty In the last digits. 
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Band 

Isotope 

Transition 

Band Center 

B' 

B" 

D' x 10 7 

D" x 10 7 


i 

I 

16 o 12 c 18 o 

(00°2) I + 00°0 

4639.5081(8) 

0.3623911(27) 

0.3681879(17) 

1.146 (36) 

1.170 (29) 

U 

4639.5082(3) 

0.3623894(14) 

0.36818571 

1.1693 (93) 

1.19030 

F 

4639.502 

0.362372 

0.3681771 

1.11 

1.17 

R 

4639.512 

0.362384 

0.368186 

1.18 

1.19 

C 

4639.5031 

0.36239118 

0.3681834 

1.1967 

1.1800 

M 

II 

16 o 12 c 16 o 

(31 1 0) III 00°0 

4591.1192(28) 

0.38996(22) 

0.39025(26) 

1.52 (14) 

1.30 (72) 

U 

4591.1195(28) 

0.389939(28) 

0.39021894 

1.585 (99) 

1.33281 

F 

4591.118 

0.389925 

0.39021894 

1.53 

1.33373 

R 

4591.027 

0.389927 

0.390221 

1.54 

1.33 

C 

4591.1188 

0.38992250 

0.3902215 

1.5284 

1.3373 

M 

III 

c 

1G 0 12 C 18 0 

(01 1 2) I + (01 1 0) I 

4614.7792(60) 

0.362857(54) 

0.3685824 

1.27 (60) 

1.14 

F 

4614.780 

0.362848 

0.3685824 

1.06 

1.14 

R 


0.362849 

0.368604 

1.19 

1.20 

C 

4614.7794 

0.3628716 

0.3685996 

1.2561 

1.1302- 

M 

m d 

16 0 12 C 18 0 

( 01 1 2 ) (01 l 0) T 

I 1 1 


0.363414(36) 

0.3691366 

1.41 (48) 

1.14 

F 

4614.780 

0.363356 

0.3691366 

1.06 

1.14 

R 

4614.787 

0.363376 

0.369156 

1.20 

1.21 

C 

f- 

4614.7808 

0.36339483 

0.3691715 

1.2033 

1.2559 

i ' • 

M 

: ■ 
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Table III 


Number of P and R lines which can be distinguished 
above the noise level and analyzed for each band. 


Band 

Number of P lines 

Number of R lines 

I 

57 

49 

II 

26 

30 

III c 

34 

34 

III d 

34 

34 
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Table IV (a) 

Rotationless Dipole Moment Matrix Elements 


Band 

» 

o 

X 

10 5 Debye 


- 

Spectrum I 

Spectrum II 

I 

136 (1) 

137 (1) 

II 

1.21 (9) 

1.19 (3) 

III 

147 (4) 

149 (4) 


Table IV (b) 

Band Intensities compared with the AFGL listing (Ref. 7) 


Band 

S° x 

22 -1 

10 (cm /mol -cm 

-*> 

i 

Spectrum I 

1 

1 

Spectrum II 

AFGL 

16 o 12 c 18 o 

0 0 
00 2 «■ 00 0 

0.122 (1) 

0.124 (2) 

0.1302 

16 o 12 c 16 o 

(31 1 0) IIJ «- 00°0 

0.019 (6) 
l = 0.164 (4) 

0.019 (2) 
5 = 0.168 (2) 

0.00205 
% = 0.18 

16 o 12 c 18 o 
(01 1 2) I «- (01 1 0) I 

0.0108 (5) 

0.0111 (6) 

0.01042 
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Table V 


Self -broadened Lorentz width parameters 
for one atmosphere, 296 K. 


1 

Parameter 

W i 

Spectrum 

Hoke and Shaw (Ref. 18) 

I 

II 

2614 cm -1 

2757 cm -1 

2 -1 
W x 10 cm 
o 

8.2 (2) 

8.34 (9) 

8.34 (2) 

8.0732 (8) 

x 10 2 cm -1 

4.8 (1.0) 

6.8 (2) 

5.83 (6) 

6.15 (4) 

W 2 

-0.15 (4) 

-0.171 (9) 

-0.1539 (4) 

-0.179 (9) 

3 -1 





V 3 x 10 cm 

3.3 (3) 

3.1 (1) 

1.82 (2) 

1.89 (9) 
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Table VI 


Values of maximum optical path difference (retardation) 
obtained by various analyses. 


REF 


Hoke (Ref. 4) 
This work 


Retardation (cm) 

Spectrum 

I 

II 

16.51 (10) 


16.64 (10) 

16.48 (15) 


16.47 (9) 
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Figure Captions 
Fig. 1 

-1 

The portion of Spectrum I analyzed. Lines of the bands centered at 4639.5 cm 
-1 

and 4591.1 cm are easily identified but lines of the much weaker band, 

-1 

centered at 4614.8 cm , can hardly be distinguished above the noise. 

Fig 2. 

Corresponding portions of the residual spectrum and the noise spectrum are 
shown. The upper spectrum is the part of the experimental spectrum shifted by 
0.2 units with respect to the calculated spectrum. 

Fig. 3 

Comparison of line position differences. Differences are plotted against 
rotational index m. The letters labeling the plots are the same as those used 
in Table II. 


Fig. 4 


1 


Calculated Intensities of the lines of band (31 0) 
plotted against rotational index m. 


Ill 


0 

«- 00 0 of 


16 12 16 
0 C 0 


Fig. 5 

The |m| dependence of the line widths for Spectrum I and Spectrum II. All 
widths have been reduced to 296 K and 760 torr. 


Fig 6 

Comparison of the |m| dependence of the widths of lines of Spectrum II with 
other results from (Ref. 18). All widths reduced to 296 K and 760 torr. 
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